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Abstract 

We present a detailed evaluation of the total semileptonic B meson width in terms of \Vcb\ 
and heavy quark parameters (quark masses and the expectation values of local heavy 
quark operators). Special attention is given to perturbative corrections which can pre- 
cisely be calculated in a scheme with a hard Wilsonian cutoff at a scale around 1 GeV 
appropriate for the OPE, and to the potential impact of higher-order power corrections. 
We point out that the latter require control over possible contributions from four-quark 
operators containing charm quark fields. Analytical expressions are given which allow 
evaluating the width with various choices of parameters; ready-to-use expressions show- 
ing the dependence on the heavy quark parameters are presented as well. We illustrate 
these results by commenting on how these parameters can be extracted and what accuracy 
is likely to be achievable in the near future. 
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Extracting the value of \ Vcb\ from the observed semileptonic width of B mesons constitutes a method 
of impeccable theoretical pedigree [1], although it often has been imprecated in the past. The theoretical 
part of the analysis proceeds in two steps: 

1. Through a heavy quark expansion (HQE) one expresses r{B IvXc) as |V^6p times a function 
of the heavy quark (HQ) parameters like quark masses and hadronic expectation values, and the 
perturbative contributions. 

2. One determines the values of these HQ parameters from observables other than the total semilep- 
tonic width. 

The first step has now reached a mature stage after having been subjected to considerable scrutiny for 
several years. The primary motivation for this note is to present an updated description of this first 
task. We will give analytical expressions as much as reasonably possible, state explicitly the proper field 
theoretical definitions of the HQ parameters and address the potential sources for theoretical uncertainties 
that arise in such an expansion. We aim at giving an 'open source code' which enables the dedicated 
reader to check the validity of the theoretical results and insert personally preferred values for the HQ 
parameters and incorporate additional constraints. 

To extract a value \Vcb\ from the observed width one has, of course to complete step two as well. 
There exists considerable information and constraints on these HQ parameters. We will sketch how they 
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can be determined. However, more work is still needed. We will present numbers in this context; yet 
they are to be understood as illustrative rather than final. 

There is already highly nontrivial evidence that the systcmatics are indeed under control. We an- 
ticipate that in the end one will be able to reduce the overall theoretical uncertainty to the one percent 
level. This benchmark will guide us in calculating rsi(i3). 

The remainder of this note will be organized as follows: in Sect. 1 we provide a summary of the 
expressions that one can use to obtain the value of \Vcb\ from the total semileptonic B width. In Sect. 2 
we present the theoretical expressions used to derive the total semileptonic B width and discuss all 
ingredients in the corresponding subsections. Wc include the effects of fom-quark operators with charm 
field previously ignored, and analyze the potential effect of higher order power corrections. Then we 
discuss in more detail the definitions and basic properties of the heavy quark operators, and address 
the perturbative corrections in Sect. 4. Section 5 summarizes the theoretical uncertainties in evaluating 
Tsi- In Sect. 6 we illustrate how one can independently determine the values of the required heavy 
quark parameters from experiment. Technical details are relegated to Appendices, where we give explicit 
expressions for the normalization point dependence of the heavy quark parameters, provide formulae 
for computing perturbative corrections implementing Wilsonian OPE, and briefiy describe the BLM 
summation to arbitrary order within this framework. 



1 Executive Summary 



We base our numerical analysis on the following expression for the semileptonic B width [2] through 
order I/toq: 



r,i(6^c) 



2(I-r)4 



zo(r) 1- 



2m2(/x) 



+ dir] 



+ ... 



(1) 



where Zo{r) is the tree-level phase space factor and r = ml{fj,)/ml{n): 

ZQ(r) = 1 - 8r + 8r^ - - 12r^ In r , 

while the expression for d{r) follows from Ref. [3]: 

34 32 32 10 

d(r) = 81nr+- - 8r^ + — - — ~ -18.3 at v^=0.25 . 



(2) 



(3) 



The electroweak correction, Aew corresponding to the ultraviolet renormalization of the 
Fermi interaction is well-known [4]: 



1 + ^ 1 + - In 



a , Mz 



TT rrib 



1.014. 



(4) 



An auxiliary scale /i is introduced to demark the border between long- and short-distance 
dynamics in the OPE. Unless stated otherwise, we adopt // ~ 1 GeV. 
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The quantity A^'^^^ accounts for the perturbative contributions; it has been calculated 
to all orders in BLM corrections and to second order in non-BLM corrections as discussed 
in detail in Sect. 2. For ^^=0.25 and = 1 GeV we find that A^^'^ ~ 0.908. 

The quantities //^, /i^, pj^ and p^g denote the expectation values of the kinetic, chro- 
momagnctic, Darwin and spin-orbit operators, respectively. Only p|, has a noticeable 
impact at 0{l/ml). In that order in l/rrii, there arise also contributions corresponding 
to four-quark operators of the generic form (6rc)(cr6). Their nonperturbative B expec- 
tation value, while small, will not vanish. For now we include such a contribution in our 
estimate of the theoretical uncertainty, given below. 

Using the numerical values mb{l GeV) = 4.6 GeV, fil{l GeV) = 0.4 GcV^ GeV) = 
0.35 GcV^ pf)(lGeV) = 0.2 GeV^ and pi5(l GeV) = -0.15 GeV^ , for = 0.25 and 
as{mb) = 0.22 we find the following expression for \Vcb\'^ 



1 , 1 



(if^)'(l-4.8[Br(i?^X„M -0.0018]) (1 + 5,,), (5) 
with the following dependence of \Vcb\ on the various heavy quark parameters: 
^^"''^ = (1 + (5th) [1 + 0.30 (a,(m6)-0.22)] x 



0.0417 

[1-0.66 (m6(l GeV) -4.6 GeV) + 0.39 (mc(l GeV) -1.15 GeV) 
+ 0.013 (/x^-0.4GeV^) + 0.09 (/9l)-0.2GeV^) 

+ 0.05 (/i^-0.35 GeV^) - 0.01 (pi^+O.lS GeV^)] . (6) 



Beyond the uncertainties in the numerical values of these parameters there is the 
potential error in |V^b| due to the limited accuracy of the theoretical expression for the 
semileptonic width we relied upon; this is denoted by 5th, for which we estimate 

5th = ±0.005pert ± 0.012hwc ± 0.004hpc ± 0.007ic . (7) 

The terms here represent the remaining uncertainty in the Wilson coefficient of the lead- 
ing operator bb (perturbative correction), as yet uncalculated perturbative corrections 
to the Wilson coefficients of chromomagnetic and Darwin operators, higher-order power 
corrections including violation of local duality in rsi(6— >c), and possible nonperturbative 
effects in the operators with charm fields, respectively. 

While the numbers of Eq. (7) are small, we consider them rather on the conservative 
side. Yet they should be viewed as what they are, namely estimates based on best efforts 
rather than mathematical theorems. 

^The number slightly differs from Ref. [5] since we adopt the Wilsonian renormalization convention 
for Darwin term as well; our p'jjin) should be distinguished from of Refs. [5, 6]. 
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2 r(B^Xc^i^); Master Formulae 

The operator product expansion (OPE) yields inclusive heavy quark decay rates as an 
asymptotic series in inverse powers of the heavy quark mass. To be more precise, there 
are different mass scales in the problem; the most important one is the energy release 
given by Er = mb — rric ior b ^ c [1]. Hence, the width in the OPE is an expansion in 
inverse powers of Ej..'^ It is useful to note that for actual quark masses the dependence 
on m.}, and rric is well approximated by rsi(-B) oc ml {mb—mc)^. 
Through order l/mg one has the general expression 



zo{r)[l+A^\r;^^)] 



\ 



2mi{ii) 



+ 327r2(l + ^Pf (r))(l-V^)2-^ + 327r2iPf (r)(l-V?)2 




(8) 



The phase space factor zo{r) has been defined in Eq. (2) and the tree coefficient for the 
Darwin term d{r) in Eq. (3). A number of clarifications are to be made about the above 
master expression. 

Following Wilson's prescription for the OPE an auxiliary scale has been introduced 
to separate the effects from long- and short-distance dynamics. Since observables do 
not depend on it, quark masses, radiative corrections and hadronic expectation values 
(S I Oils) combine to yield a compensating /i dependence as indicated in Eq. (8). The 
natural and most efficient choice is // ~ 1 GeV, which is particularly relevant for the quark 
masses [8]. 

It should be kept in mind that the B state over which the expectation values are 
taken is the real B meson, not the meson in the heavy quark limit m;, — > 00. In this 

convention the total width to order is corrected by a single Darwin expectation 

value (modulo the effect of four-quark operators to be addressed below). No non-local 
correlators enter, and the spin-orbital expectation value p|g enters only as a l/m^ piece 
of the total Lorentz-scalar chromomagnetic expectation value [9].^ Likewise, the kinetic 

^While the OPE in Eqs. (1), (8) at first siglit look like expansion in powers of l/mi,, this is not really 
the case. In the SV limit mb~mc <^ rUb the expansion parameter is obviously l/{mb — Wc); at m.b^m.c 
the expansion becomes meaningless no matter how large 777,5 is- Constructing the OPE directly for the 
widths shows that the expansion parameter in general is 1/ (mb — nic) [7], which comes about in Eqs.(l), 
(8) due to the interplay of l/rUb and r = m\/m\. On the other hand expectation values like {B\pb\B) are 
expanded in powers of l/mb- 

■'The inclusive moments of various decay distributions do exhibit explicit dependence on p\g] it is 
however too weak to matter in practice. 
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expectation value enters solely through the l/m^ expansion of the leading operator bb 
[10, 2, 9]. 

All Wilson coefficients are given by short-distance dynamics and in practice are evalu- 
ated in perturbation theory. In the convention adopted in Eq. (8) the tree-level coefficients 
zq, d etc. are modified by perturbative corrections A^^^^ starting with at least one power 
of ag. 

The leading nonperturbative corrections arise in order l/rriQ and are controlled by the 
expectation values and fJ^cif^) of the kinetic and chromomagnetic dimension- five 

operators {Di, denotes the covariant derivative), respectively 

(Since wc have explicitly introduced the effect of the LS operator, fiQ above should be 
understood as the average of only the chromomagnetic field, {(^B), although it includes 
l/nib finite-mass effects.) The Darwin and spin-orbital LS terms p^^iJi) and PLsif^) emerge 
from dimension-six operators: 

'"^^'"^ ^ ^^^'^^-^^■^^^l^^'^ ' P'^^^l"^ ^ ^{B\b{a-ExzD)b\B), . (10) 

The Wilson coefficients for the operators describing nonperturbative effects are completely 
known only at tree level. Since they represent power-suppressed effects showing up at the 
one percent level, this is a reasonable approximation. However, it is desirable to improve 
it in the future computing 0{as) corrections to them. 

The last term in Eq. (8) proportional to Fg denotes the effect of generic SU (3)-singlet 
four-quark operators, other than Darwin operator, of the form hVbqVq with the sum over 
q = u,d, s, and T including both color and Lorentz matrices (to the leading order in l/rrib 
only 7o x or 7^75 x 7*75 structures survive, but one does not need to rely on this). 
Their Wilson coefficients are order 0{as), and we neglect these contributions. 

On the other hand, we have explicitly included in the master equation a possible effect 
of the (tree level) expectation values of the four-quark operators with the charm field. It 
is routinely skipped assuming that it vanishes due to the sizeable charm mass. While 
it represents a reasonable approximation, it can be valid only up to a certain accuracy. 
Since the Wilson coefficient of one linear combination of four (in the heavy-6 limit) such 
operators emerges at tree level and is strongly enhanced by the effectively two-body phase 
space, this effect might not be totally negligible. To complement its purely theoretical 
estimates or upper bounds, we find it advantageous to explicitly introduce this effect and 
study its possible manifestations in experiment: 

He = ^{B\hu{i-i,)cmi-i5)m, Hp+^'^^p) , = ^ ■ (11) 

This way to refer to such nonperturbative effects literally applies when normalization 
point n is taken above rric. Having the normalization scale below rric amounts to the c 
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quark field being integrated out; then there are no dynamical c quark fields. Nevertheless, 
the physical efi^ects lumped into still do not have to vanish exactly; they are rather 
partially reshuffled into higher-dimension b-quark operators with only light fields. 

The second charm-related D = 3 parameter He in Eq. (8) represents the cumulative 
effect of the remaining combinations of the four-fermion operators with charm quarks. It 
is expected to play an insignificant role, but will be briefiy addressed as well. 

It turns out that the consistent 1/mf, expansion to higher orders mandates introducing 
such four-quark charmed operators. They can be excluded only at the price of spoiling 
the expansion: instead of l/m?, or 1/Ej. some of the higher-dimensional operators would 
get powers of rric in the denominator to compensate for the increasing dimension. To 
state it differently: the Wilson coefficients of the higher-dimensional operator can contain 
an enhancement by positive powers of 1/r. The aggregate effect of such terms can be 
viewed as the definition of He at low normalization scale. Moreover, a significant fraction 
of the numerically enhanced Darwin coefficient function d{r) represents actually the mid- 
virtuality (semi-hard, or hybrid) contribution to the leading charm operator. A detailed 
discussion of these operators and their relation to the notion of "Intrinsic Charm" will be 
given in a separate publication [11]. We shall briefiy recapitulate the main points of that 
analysis below in Sect. 3.3. 

The Wilsonian separation of "soft" and "hard" effects can be done in different ways. 
In particular, this affects the precise values of the heavy quark masses at a given /i. We 
consistently use the scheme based on the Small Velocity, or Shifman-Voloshin (SV) sum 
rules which introduce the normalization scale ^ via the cutoff over the excitation energy 
of the hadronic states [8, 12]. Other possibilities to define heavy quark masses discussed 
in the hterature are considered in the review [13].^ The advantages of this choice are 
transparent. First, they allow to define the heavy quark masses as well as the expectation 
values of the higher-dimensional operators on a parallel footing. Secondly, being expressed 
via (in principle) measurable quantities, they have a definite value in QCD at any given 
normalization point. The question of their precise extraction from experiment then be- 
comes meaningful (see, e.g. Ref. [14]), without invoking doubtful procedures. Thirdly, 
such a renormalization scheme respects exact inequalities established in the heavy quark 
limit (for a recent review, see Ref. [7]), which turn out very constraining in practice and 
allow to go far beyond naive dimensional estimates even for higher- dimensional operators. 
Finally, it has been demonstrated in numerous applications that using such a scheme typ- 
ically yields a much improved convergence of the perturbative expansion, see, e.g. a recent 
experimental implementation in Ref. [6]. 

Below we discuss the salient features in more detail. 

^Their assessment in Ref. [13] is not always justified in our opinion. 
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3 Heavy Quark Parameters 



3.1 Heavy quark masses 

With the width depending on the fifth power of the heavy quark mass any uncertainty 
associated with it - conceptual or numerical - has a grave impact on the accuracy with 
which \Vcb\ can be extracted. Therefore great care has to be applied to properly defining 
it and consistently calculating the radiative corrections. 

The pole mass is defined in analogy to the electron case in QED as the position of the 
pole in the quark Green function. This pole mass is gauge invariant and perturbatively 
infrared finite. Employing it is often convenient for purely perturbative calculations. 
Yet in full QCD due to confinement there is no such pole - the pole mass is thus not 
'infrared stable' due to nonpcrturbativc dynamics. Moreover, it has an intrinsic and thus 
irreducible theoretical uncertainty ~Aqcd already in perturbation theory; this makes it 
inappropriate when one wants to include nonperturbative contributions, which are power 
suppressed [15]: 

r,i(S) oc ~ {mfy{l + 50(AQCD)/m6) ; (12) 

i.e., the intrinsic uncertainty is at least parametrically larger than the leading nonpertur- 
bative contributions that arise in order 1/m^. The meaning of this observation can be 
easily understood. A long-distance pole mass in perturbation theory is a counterpart of 
the heavy flavor hadron mass of full QCD, since it includes effects of gluons with arbi- 
trarily small momenta. In this respect the analogue of the quark mass in QCD where 
nonperturbative effects exit, is a short-distance mass mq^ix) of perturbation theory, since 
such a mass excludes soft gluons interacting strongly. The central OPE result that there 
are no contributions ~ C(l/m6) contributions to fully inclusive widths [10, 16] - in other 
words, the hadron {B meson) mass is irrelevant for the width which depends only on mf,, 
but not on Mb— mi, - translates then into the statement that short-distance masses have 
to be used when aiming at the power-like accuracy. 

To a limited extent, the pole mass can still be used to approximate the width: the 
order-to-order shift in the apparent position of the pole quark mass in perturbation theory 
is offset by a significant change in the overall perturbative correction factor. Yet such a 
procedure may not allow to properly account already for the actual leading nonpertur- 
bative corrections. This is illustrated by the following simple observation. Starting with 
the short-distance mass, one has for the leading (parton) contribution 

rsi(5)ocAP^'*(/x)m^(/i) (13) 

with neither factor having uncontrollable infrared pieces. The same perturbative expres- 
sion for the pole masses then must be of the form 



OC AP-*(;.)(1-5A) (^Pole)5 ^ ^P-t ^5(^) 



2 

1 - 15 ■ ' ' 



(14) 



where A/mj expresses the uncontrollable infrared contributions in the ratio of the pole 
to short-distance mass. The contribution ~ C(l/mb) has indeed disappeared, yet the 



7 



expression is still deficient: while there are contributions ~(9(l/m^), they have to be of 
the form 1 — regardless of the particular interactions.^ Hence the last term in Eq. (14) 
coming with the coefficient 15 has no counterpart in actual QCD and is rather an artefact 
of using the pole masses. 

An example of a short-distance b mass is the MS mass inb{mb). The normalization 
scale m,ij however is unnaturally high for the problem; therefore it introduces significant 
higher-order corrections not related to running of ag per se, which are under poor control 
and thus limit the accuracy [8]. 

A better choice when treating B decays is taking the Wilsonian factorization scale /i 
around 1 GeV. An example of such a mass is the one defined through the SV sum rules, 
the prescription valid to arbitrary perturbative order. It has the meaning of the mass 
entering the kinetic energy of the heavy quark and is often called kinetic mass m^^^{fi). 
At every given it represents an observable physical quantity and, with an appropriate 
choice of fj, its extraction from experiment provides good stability with respect to radiative 
corrections. Its definition is discussed in Sect. A.l. Here we mention that its value has 
been extracted from the threshold domain of the e^e~ annihilation into bb , 

m^^"(l GeV) ~ (4.57 ± 0.06) GeV (15) 

Prom now on we omit the superscript "kin" using this mass as the default definition. 

For the charm quark, the separation scale taken around 1 GeV is similar in magnitude 
to rrif. itself. The kinetic mass mf.{ij,) can still be used, with the realization, however, that 
this definition starts loosing direct physical meaning for increasing ^. 

Historically heavy use has been made of relating the quark mass difference nib — mc 
to the spin-averaged meson masses to reduce the number of parameters and to suppress 
ambiguities related to heavy mass definitions: 

rrib-m,^ Mb -Md + ^ (—-—) + f A " + Oil/m'g) (16) 

2 Kmc nibJ 4 \mi ml) ^ ' v 

where Mb[d\ = Mb[d\/^ + 3MB*p*j/4 denotes the spin averaged meson masses and p^ — 
+ pI the sum of two positive (in the advocated scheme) nonlocal correlators (they are 
defined in Ref. [9]). Eq. (16) has become a widely used relation: for once one realizes that 
p"^ cannot differ from 0.45 GeV^ by significantly more than 0.1 GeV^, it leads to a reduced 
uncertainty in m^— mc. Having the same scheme and normalization point for both mij and 
mc has an advantage of making the heavy quark expansion relation between mb—mc and 
Mb — Md more direct. If such constraints are not imposed, it is justified to use ifidmc) 
as well. We argue in Sect. 6 that it is safer not to invoke these relations. 

3.2 Operators of dimension 5 and 6 

The precise definition of higher-dimensional operators affects perturbative corrections to 
a lesser extent than the choice of the heavy quark mass scheme. Yet it is required to 

^The perturbative contributions discussed here cannot give rise to spin-dependent effects contained in 
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make sense of assigning them a definite numerical value; once again the values have to 
be normalization-scale dependent. In particular this refers to the spin-singlet expectation 
values and pjj. They can be defined in the same physical way as the heavy quark 
masses; the formal definition is given, e.g. in [14]. Once given, it allows to sensibly 
extract their values; for instance, Ref. [14] concluded from the hyperfine mass splitting 
that 

fi^il GeV) = (0.35^:°^) GeV^ . (17) 

Kinetic expectation value /j,^ a priori is less certain, however in this regularization scheme 
the inequality 

f^'M > i^Ui^) (18) 

holds for any normalization scale. With a number of experimental data strongly favoring 
GeV) ^ 0.5 GeV^ a limited range is left allowed, and the estimate 

I^Kl GeV)) = (0.45 ± 0.1) GeV^ (19) 

seems rather conservative. It is also supported by purely theoretical estimates including 
QCD sum rules [17]. ^ 

A less precise estimate can be obtained for the Darwin expectation value, p%{l GeV) ~ 
(0.2 ± 0.1) GeV^, yet it must be positive. The spin-orbit average p^^g is expected to be 
negative. It obeys the constraints 

-pIs < pI, \pIs\ < 2PD ■ (20) 

Since the effect of this operator is typically strongly suppressed compared to the Darwin 
term, already these general bounds suffice for practical purposes. 

A clarification is in order here. Strictly speaking, all the above bounds hold in the 
heavy quark limit, i.e. when l/rrib corrections in the expectation values arc absent. It 
has been argued [14] that their shift from the nib oo limit is negligible in B meson. 
Not only it is governed by the parameter l/2m6, it should be additionally suppressed by 
virtue of the small excess of p^ over Pq manifesting proximity to the so-called BPS limit 
for the heavy meson ground state. Then the deviation of these hadronic parameters lies 
well below available precision, and can be neglected in practice. 

3.3 Operators with charm quarks 

As mentioned in the Executive Summary and the beginning of Sect. 2, four-quark op- 
erators containing a pair of charm and anti-charm quark fields necessarily arise in the 
OPE and their B meson expectation values do not vanish exactly. These effects will be 
discussed in a forthcoming dedicated publication [11]. Here we give a brief introduction 
and state the salient conclusions of that analysis. 

In discussing such effects in the OPE framework, it is essential to distinguish between 
an expansion in I /nib (actually, l/{nib—nic)) and in l/nic- For the sake of clarity we resort 

^According to the FNAL group, lattice studies suggest so far a somewhat larger value [18]. 
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to a somewhat idealized scenario where the hierarchy between the two scales is stretched 
up compared to actual QCD. 

To the leading third order in l/m^ there are four such operators of which two are 
relevant for rsi(6— >c): 

^'^^^^^^ ^ -^{B\H{l-l,)ccf{l-^,)b\B) 

Fc = m - -^(5|^f%(l-75)ccfy(l-75)6|i?). (21) 

(Fc is the tree-level counterpart of He in Eq. (8)). To higher orders in l/m^ there are 
additional operators with derivatives acting on the quark fields. They can be reasonably 
neglected in the total width. For general purposes it is convenient to consider the extended 
set of four such operators 

Oy — bb CTfQC 0\ = bab C775C 

= bfb c^joc 0\ = b^ab c^'fj^c (22) 

which define the expectation values He and Fc by virtue of Fierz transformations: 




The indices V, A and s, o correspond to products of vector or axial, and color-singlet or 
-octet flavor-diagonal bb and cc currents. Below we shall mainly refer to the contributions 
related to H^. as most relevant for rsi{b—>-c). 

The existence of such an effect can be understood since this expectation value simply 
states that the propagation of the decay quark inside B meson is not totally perturbative. 
It is related to local expectation values because we consider the inclusive width coming 
from c quark momenta much smaller than the energy release where the momentum of the 
lepton pair does not differs much from = {Ej., 0, 0, 0), and we expand in this difference 
[19]. This inclusive probability for transitions to charm quarks is affected by nonpertur- 
bative dynamics to a lesser degree than for decays to light quarks; yet it is evident that it 
should be present at some level. Its magnitude can be qualified and estimated [11]. The 
relation of these effects to the "Intrinsic Charm" in B mesons is discussed there. 

The charm quark propagator undergoes nonperturbative effects even when it is hard or 
has large energy or momentum. These effects, however are accounted for in the standard 
OPE; to the leading order they are given in terms of /x^ and Hq as exemplified by the 
sum rules considered in Ref. [9]. The effect of the Darwin term comes from both hard 
and soft charm configurations. Hence, in the Wilsonian approach a part of the effect of 
the Darwin operator in Eq. (8) actually represents the contribution H^ of the four-quark 
charm operator. This part belonging to p% depends on the choice of the normalization 
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point yU. and vanishes if jj, is taken well below tUc - this option is possible only for a heavy 
quark. 

If the charm field can be considered as heavy in the scale of typical hadronic masses 
/^hadr, One Can apply heavy quark expansion to these nonperturbative expectation values. 
It proceeds, however, in 1/mc and not in l/mj,: 



Similar to the OPE for inclusive widths, it includes only local h...h operators which, 
however, now involve only light (gluon and rt, d, s quark) fields. In the above equation 
denotes the dimension of the light field operator Ok- Therefore, without introducing the 
explicit operators with charm, there exists a chain of higher-order power terms scaling 

A" A" 

like o*^*?.. rather than . It is instructive to keep in mind that the large coefficient 

of the Darwin operator is associated with the lowest term in this subseries corresponding 
to n = 3 which gives rise to the In — enhancement. 

The l/rric expansion of the charmed expectation values thus requires classifying pos- 
sible operators Ok and computing the coefficients Cjk-'^ It is shown [11] that at the tree 
level there arc no contributions to order 1/mc, while l/m^ effects are driven by a few 
operators, depending on the particular color and Lorentz structure.^ The dominant in Nc 
structure for semileptonic decay is color octet, see Eq. (23), for which l/m^ efi'ects are 
present for both vector and axial currents. 

It is interesting, however, that absence of l/rric nonperturbative charm loop effects 
can be vitiated once pcrturbativc corrections arc included [11]. There are two spin-singlet 
operators Oe'^ and Oe-e and one spin-dependent OexE operator which emerge already to 
order ^^(mc) from the dominant (in semileptonic decays) color-octet operators, vector 
and axial, respectively: 

OEE = ^^h{l-^Ti)E.E 6(0) = ^ 6 (1 - ^Tr ) (Tr.TroTroTrfc) h 

Oexe^,^ iba-ExEb{0) = -^ie''"''^ bj^j^n^TTpTTonxb , 

where E is the non-Abelian chromoelectric field strength (it includes coupling gs), and tt^ 
are the nonrelativistic energy-momentum operators. An example of the Lorentz-invariant 
operator yielding such spin-dependent interaction is 



^Alternatively, one can integrate the sum rules for the corresponding combinations of the zeroth 
moments (for the leading-mfe operators) of the heavy quark structure functions studied in [9], over q. 
This would require, however extending the expansion at least to two more orders in l/rric- 

^Say, no contribution for the color-straight operator with vector currents, one for the product of 
color-straight axial currents, etc. 
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with denoting the full-QCD covariant derivatives. Such operators come with coefficient 
oc as{mc)/7i and are not expected to dominate for actual charm quark. They would 
represent the leading effect in the true heavy charm limit. Physically l/rric corrections can 
describe interference between the (gluon-dressed) parton decay amplitude and internal cc 
component of the meson, see Fig. 1 similar to interference discussed in detail in Ref . [20] . 
Soft fields shape the internal charm component of wavefunction, while short-distance 
effects here are characterized by virtuality starting at rric and governed by as{mc)- To 
even higher order in as{mc), and/or including gluons with momentum scaling as more 
operators become possible. 



1/ 




q q 

Figure 1: Example of the interference contribution to the decay width suppressed by a single 
power of l/rric- Shown gluon is hard. All cuts should be included to maintain the correct scaling. 

The effects of intrinsic charm in the context of DIS on nucleons were thoroughly inves- 
tigated in Ref. [21], though without perturbative effects. The relevant charm operators 
are quark bilinears in this case, yielding quite different pattern of the power expansion. 
Including hard perturbative corrections can change leading power in l/rric for heavy quark 
decays, as noted above. 

The practical OPE Eq. (24) in the form of an expansion in Aqcd/'^c is only asymp- 
totic; the charm expectation values contain also "exponential" pieces like e~^™"=/^''="*'. 
Various considerations lead us to expect that for actual charm with mc — 1.2 GeV these 
exponential terms can be, if not dominant, at least not too far suppressed compared 
to the leading power like effects. These contributions would represent the clearest and 
unambiguous example of genuinely independent effects associated with nonperturbative 
corrections for charm quarks in B mesons. 

All these contributions, which are suppressed by powers of l/rric and therefore could 
be sizeable thus posing a danger to the l/mq expansion, just reflect the existence of 
operators in the OPE of the width containing pairs of c and c flelds, namely Oj (and 
their higher-dimensional analogies). Their impact is best analyzed directly through the B 
expectation values of Oj. The size of these nonperturbative expectation values is however 
not very certain. For rsi(6— >c) there is a single relevant operator, Eq. (11). The size of 
its expectation value \Hc\ is estimated [11] 

[i/cl ^ 0.005 GeV^ (26) 
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which translates into a possible contribution to the width up to 1 to 1.5%. This looks 
reasonable if we compare it with the effect of the Darwin expectation value itself which, 
according to our master equations is estimated to be about 3%. A fraction of the Darwin 
term effect, on the other hand, represents the leading, 1/m^ piece of the four-quark 
operators with charm. 

An alternative way to arrive at a similar numerical estimate is calibrating them with 
the size of nonfactorizable contribution of four-quark operators bbqq with light quarks q, 
which would enter the total semileptonic b ^ u width. The four-quark operators with 
non- valence hght quark were estimated in Ref. [22] to affect the width at a couple percent 
level, employing the information about D meson widths. Assuming that the expectation 
values with the charm field replacing light non-valence quarks exhibit extra suppression 
due to the charm mass by a factor of four, and accounting for the difference in the phase 
space in the two decays, we end up with an estimate somewhat below the one percent 
level in r^i{B). 

3.3.1 Experimental constraints 

A more robust approach is to constrain the charm-related expectation values directly 
from experiment. Available studies of the semileptonic distributions already yield some 
more indirect bounds on such effects, however the limits fall quite short of the scale 
expected theoretically. One reason behind this insensitivity is transparent: they manifest 
themselves in the same way as the soft part of the Darwin term, and therefore to some 
extent can be absorbed by redefining the effective value of p^. In particular, for sufficiently 
inclusive characteristics replacing 

dpl^dpl + 32Tr\l-V^YH, (27) 

may account for the bulk of the effect, and the actual difference in the effect of Darwin 
operator and the charm-related nonpcrturbativc four-quark expectation value will be 
revealed at a suppressed level only in more subtle quantities, say in the dependence of the 
hadronic distributions on lepton energy. 

There can be various physical manifestations of the nonperturbative effects associated 
with moderately short-distance charm quarks in non-charm beauty hadrons. A number 
of effects have been discussed in connection to the "intrinsic charm" hypothesis [23]; it is 
not always easy to disentangle them from other underlying effects, however. 

Here we are concerned with a more concrete aspect, namely the way such effects may 
manifest themselves in the inclusive semileptonic decays via local 6-quark operators with 
charm, Eqs. (22). The search for an unambiguous manifestation is complicated by the 
expected percent scale in rsi(6 c). Yet even any higher direct upper bound would 
constitute a very valuable information complementary to theoretical estimates. 

In the heavy quark expansion these four-quark operators would most directly show up 
as 5-like contributions near the quasi-two-body kinematics b — > slowc + iu, which naively 
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corresponds to such leading terms in the OPE, with 



1 dr,, _ igtt^ h 

6 dq^-{MB-MD*) 



Tsi 


dg2 
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dr^c 




dEx 
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dVcc 




dEe 



3 



Mb-Md* ml 

• 327r^ ^S(Ex-{Mb-Md*)) , 

mfj 

: 3271^^6 (Ee-M^^) , (28) 

mf ^ '' 

etc. Placed in the narrow slice of the whole kinematic domain, these effects are locally 
enhanced and potentially detectable in differential distributions, much in the same way 
as was suggested for nonfactorizable effects in bbqq-type operators with light quarks in 
the lepton spectra [19], or Mx or Ex distributions [24]. 




Q Q 

Figure 2: Nonperturbative charm-field effects in the OPE for the inclusive decays. With charm 



fields being low-momentum, the kinematics is close to two-body for leptons, with c^mij—mc ■ 



There are also similar operators with derivatives. In the total widths they are sup- 
pressed by additional powers of l/m^. They lead to smearing the above quasi-two-body 
peaks over the domain of typical width ~ Aqcd- The smearing in practice can stretch 
over a sizable fraction of the allowed kinematic range for actual m^.^ It should be noted 
that the charm operators can lead to both enhancement or depletion of the decay rates, 
depending on the sign of the expectation values. 

The above pattern suggests the strategy of studying the integral over the particular 
low-E'x kinematic domain and comparing it with the expectations based on the parameters 
extracted from the total moments. It turns out that the overall integrated hadronic energy 
{Ex) and the invariant mass {M\) are theoretically expected to depend on practically 
one and the same combination of parameters mi) — 0.67mc + O.lyU^ — 0.25pf3. Comparing 
the two averages experimentally would then be a cross check of possible additional effects, 
including charm expectation values. 

^Pcrturbative corrections associated with the anomalous dimensions of the leading operators yield 
slower decreasing tails as discussed in Ref. [19] ; yet in practice this perturbative smearing should hardly 
be relevant. 
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It seems to us that the available accuracy in the fully integrated moments (with respect 
to Ei, Mx or q^) does not allow at present to constrain the expectation values, in particular 
He with the precision even nearly matching the theoretical expectations Eq. (26). An 
interesting possibihty is to study the partially integrated over Ex double differential 
distribution 

For the four-fermion operators it is expected to exhibit a more or less flat contribution 
above a hadronic scale Mq for moderate g^, 



j dg2 A(g2) ~ 32tt'' 



ml 



which however fades away with increase of q"^ . In this scenario the effect can originate 
from a limited range of Mx below a hadronic mass M^* + Mq with Mq expected to be 
around mc and/or //hadr ~ 

IGeV. 

The most precision data available now are from CLEO lepton moments measured 
with a lower cut 1.5 GeV [25]. However, the average lepton energy in this interval is 
{E() ~ 1.73 GeV which is very close to the expected location of the discussed contributions 

ET ~ ^ 1.64 GeV . (30) 

This strongly suppresses the sensitivity of these moments to H^.- At the same time, this 
physical example illustrates the fact that the integrated widths with cuts are far more 
vulnerable to higher-order effects: the relative sensitivity of just the rate with Ei > 
1.5 GeV to He can naturally be three times higher than of Tlf{b^c) simply because such 
a partial width comprises only about 30% of the total rate, and impact on the CLEO's 
Rq can be even more significant. Yet we think it may be premature to attribute some 
apparent discrepancy between, say, measured values of Ro and Ri to effects of charm 
operators, in view of insufficient control over usual higher-dimension operators for the 
widths with such a cut [26] . 

It is also worth emphasizing that at a high enough cut on E^ the operators describing 
the addressed nonperturbative effects cease to remain local (in space). Whether the cut 
at 1.5 GeV is safe in this respect, must yet be understood. 



3.4 Heavy quark operators of D>7 

The important question about the achievable accuracy of the OPE for rsi(6^c) is the 
magnitude of higher-order power corrections beyond l/mf terms. Their number prolifer- 
ates, and it might look hopeless to get a meaningful answer. Nevertheless, we shall argue 
that the effect of the higher-dimensional operators while applying the proper Wilsonian 
procedure should be suppressed and not exceed the percent level. One also needs to qual- 
itatively understand the hierarchy of the computed corrections through order l/m\. This 
is discussed later in this section. 
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To get the way around apparent complexity of higher-order corrections, we use a 
few basic ideas. A few different physical momentum scales enter the problem even if one 
abstracts from perturbative corrections generally messing up physics at different distances. 
Namely, a number of corrections are driven solely by the scale l/2mb - among those are 
preasymptotic corrections to the expectation values of nonpcrturbativc operators over the 
actual B meson state. It is clear that these terms among higher-order operators can be 
neglected at once. 

The second largest scale specific to the inclusive width is the energy release mb—rric ^ 
3.5 GeV. Since it still is significant, the higher-dimension operators controlled by such an 
expansion parameter yield the contribution far below the percent level as is illustrated 
by the LS term, unless their effect is particularly enhanced. There are two reasons for a 
possible enhancement of the coefficient functions. One is the "infrared instability" when 
rric/mb — > 0. It was discussed in Sect. 3.3, and they are taken care of by introducing 
the explicit expectation values He, He of the four-quark operators with charm in Eq. (8). 
Indeed, infrared sensitivity is possible only for soft c quark line, with the lepton pair 
carrying the momentum close to mb—mc- This effectively contracts the two weak vertices 
in space and time yielding local four-fermion operator bVccTb. 

The second reason discussed in detail in Ref. [8] is due to presence of the large pa- 
rameter N — 5 describing the power N with which heavy quark masses enter the total 
width, r oc {mb,mb — mcY. Therefore, to find the way through the mace of higher-order 
power corrections, we will use the classification of the coefficients over this parameter. 
Namely, we can require an extra power of the hadronic scale Hha.dr/'tTT'b in the nonpertur- 
bative corrections to be accompanied by the maximal power of N. 

The 1/N classification is quite transparent and yields a number of useful facts. As 
pointed out in Ref. [8], large N pushes one to the so-called Extended SV regime, since 
large places most of the weight on maximal of the lepton pair. It is also evident 
that this dominance is further enhanced for higher-dimensional operators: the leading 
contribution comes from multiple differentiation of the c quark propagator with respect 
to ml, which pushes it further and further into the infrared. In this respect the effect 
is the same as for the first mechanism, infrared enhancement, and in the same way the 
largest pieces are absorbed into the explicit charmed expectation values. The upshot of 
this approach is two-fold: 

• The dominant operators are spin-singlet like the kinetic and Darwin terms, as op- 
posed to the spin-dependent chromomagnetic or LS operators. These are among those 
which determine higher moments of the light-cone distribution function. The light- 
cone set of operators would dominate if rrib were really much greater than rric so that 
nib/N > rric, the case where normalization point ^ for the four-fermion operators with 
charm could be taken well above rric- 

• In the realistic case, the leading ones among these spin-singlet operators are those 
coiitaiiiiiig inoio time (ieii\'ati\'es: extra space derivatives generally yield subleading in 

^°The increase in the Wilson coefficients should be distinguished from the expected factorial growth of 
the matrix elements of the operators compared to the dimensional estimate Mhadri latter is related to 
violation of local duality [27, 28] and will be addressed separately. 



16 



effects. Since the left-most and riglit-most covariant derivatives in (bD...Db) must be 
space-like for the heavy b quark, it is actually the SV-type operators determining the SV 
heavy quark distribution function that dominate the higher-order power corrections. This 
is in a clear accord with the onset of the Extended SV regime mentioned above. 

It is worth noting that the dominance of the SV-like operators makes the estimates 
more conclusive. They can be fully defined with explicit normalization-point dependence 
in the same way as the operators through D = 6 , sec, e.g. [29, 7, 14]: 

{biD^iiDoTiD, 6), = 3E [2{e^j,ryi),\^ + ie^;^^^^-^ \'] . (31) 

e<jLi 

This eliminates the notorious problem with fixing the ad hoc factors which often can 
be arbitrarily reshuffled between the operators and their Wilson coefficients, obscuring 
numerical estimates. For instance, the parameter A2 used in some papers historically was 
defined with the factor 1/3 of the actual size of the physical chromomagnetic operator 
in B, and HQET's Ai was laxly defined as minus the expectation value of the positive 
operator. 

The chain of SV expectation values has a definite sign. In addition, their scaling is 
given by a physical mass, the energy of the excited P-wave states at which one observes 
the onset of the perturbative regime. The recent data on the Mx spectrum in semileptonic 
B decays [30, 31] suggest that this scale is near the expected 1 GeV. It is clear that the 
magnitude of the ratio of this mass to the energy release mb — mc would govern the rate 
of convergence of power expansion. For the SV regime, for example, this is explicit in the 
analysis of the Orsay group [32]. The detailed study of the Mx spectrum in semileptonic 
decays will allow even some more quantitative bounds on the leading higher-order heavy 
quark operators. 

The analysis of possible higher-order effects conducted along these lines suggests that 
the effect of higher-order operators with D>7 constitutes at most a fraction of percent 
in rsi(6^c), at least if the four-quark operators with charm field are properly introduced 
to take care of the infrared contributions. It is worth commenting in this respect on the 
numerical pattern in low-order power corrections, which is not always properly interpreted 
and may lead to unjustified suspicions about convergence. 

According to Eq. (1), the first few power corrections amount, in fractional units to 

S^2 ~ -0.012 , 6^2^ ~ -0.035 , Sp3^ ~ -0.035 , 5^3^ ~ -0.003 . (32) 

The values /i^ = 0.4 GeV^ /x^ = 0.35 GeV^ pi = 0.2 GeV^ and pis = -0.15 GeV^ have 
been used. The significant contribution of the Darwin term compared to the second-order 
corrections, and in particular, to the kinetic operator, may raise some concern. In fact, 
applying the above theoretical perspective one can realize that the Darwin operator has 
the expected, normal size. 

First of all, as detailed earlier, a significant part of the Darwin expectation value 
actually represents the contribution of the four-quark charm operator He in Eq. (11) 
if it is understood in the Wilsonian sense normalized at the appropriate scale, say at 
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jj, ~ rric + 1 GeV. The remaining part is then noticeably smaller than the estimate of 
Eq. (32). 

It is more important to appreciate the following fact: the effect of the kinetic and 
chromomagnetic operators is actually essentially suppressed compared to a priori expec- 
tations based on the scale arguments. These are nontrivial dynamic results of the OPE 
which cannot be accounted for by simple dimensional considerations. 

The most obvious example comes from l/rrib effects. Dimensional estimates would 
yield 



the result typically exhibited by naive quark model calculations. Such a correction would 
actually belong to the SV family in Eq. (31): according to Voloshin's "optical" sum rule 
I A amounts to the expectation value with n = —1. Such corrections physically exits. 
However, the OPE ensures that in QCD they exactly cancel between the binding effects 
in the initial B state and hadronization corrections in the final state. Clearly one would 
not interpret the fact that 1/m^ corrections to total widths are infinitely larger than the 
absent l/rrib effects, as a non-convergence of the power series! 

It is even more interesting that a similar cancellation, in a sense holds even to second 
order ml/ nib. Naively one would expect the kinetic operator to appear with the coefficient 
scaling like 5-4/2 ^ 1, and to have powers of m^ — mc rather than rrib in denominator. 
However, the OPE once again ensures that such an effect is totally absent from the 
inclusive width if strong interactions are described by a gauge theory like QCD.^^ The 
only effect of the kinetic operator emerges as the I /ml correction to the expectation value 
{B\bb\B) [10, 34, 9]. This makes evident why its effect is numerically suppressed. We 
can safely discard such contributions for higher-order operators, which has been already 
assumed in formulating the rules of the analysis earlier in this Section. 

The coefficient of the chromomagnetic operator docs not vanish. Yet it is still partially 
suppressed, which can be explicitly seen computing its effect when employing iV = 5 as a 
free parameter: it scales hke N, but not N'^. This is related to the spin-nonsinglet struc- 
ture of the chromomagnetic operator - it does not belong to the SV family of operators 
of Eq. (31) being an antisymmetric combination of the spacclikc derivatives. This partial 
suppression is also reflected in the dependence of its coefficient in the width: compared 
to the parton width it is enhanced by a single power of rrib/ {mb — mc). 

A similar partial suppression is present for all higher-order spin-nonsinglet operators, 
which is precisely illustrated by the LS contribution. It has a coefficient scahng like that 
of fiQ, and its numerical contribution is safely at a sub-percent level. 

For the Darwin operator, however, the cancellations 'accidentally' enforced on the 
lower-dimension SV parameters in the OPE do not hold anymore, and its coefficient 
has 'normal' or expected magnitude. Its numerical impact on the total width is on 
general grounds expected to be in a few percent range, assuming the normal scaling 

^^Even refined parton-based models typically predict a large positive effect of the kinetic expectation 
value; see, for example the light-cone based approach of Ref. [33] yielding the coefficient -1-35/6. 
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of the expectation values following from the family of the SV sum rules [1, 7]. The effect 
of higher-order operators then should not exceed the scale of a half percent. 

To summarize, we expect, and have argued that the contribution of heavy quark 
operators with D>1 can give only a sub-% effect on rsi(6^c) provided the four-fermion 
operators with charm field in the Wilsonian OPE have been properly incorporated. 

4 Perturbative Contributions 

Once the perturbative effects in the quark masses and higher-dimensional operators have 
been properly incorporated, one can turn to the radiative corrections to the Wilson coeffi- 
cients. In practice, high accuracy is required for the leading (parton) operator coefficient 
denoted as 1 + A^^^^ in Eq. (8), and we routinely refer to it as the perturbative correction 

AP^"^*, see Eq. (1). 

It has become customary to express perturbative corrections as series in the MS cou- 
pling evaluated at m^. Yet the scales typical for the radiative corrections are notably lower 
than m^. By evaluating the lowest order term in Q;^^(m(,) one generally underestimates 
^pert g^j^^ forces the a1 and higher terms to possess large coefficients. If, however, at 
least the second order BLM contributions are known, then one can express the findings 
in terms of a^^{mb) without raising numerical havoc. 

It should be noted that using the Wilsonian prescription with an appropriate hard 
factorization scale ^l significantly improves convergence of the perturbative series and 
makes the above problem less acute. This is illustrated by Fig. 3 showing the contribution 
of different momentum scales in the total width, with and without separation of short- 
and long-distance effects. The effect of removing the infrared domain is self-manifest. 
Moreover, it is evident that a too significant part of the perturbative corrections with 
pole masses comes from gluon momenta below 500 MeV, which casts serious doubts on 
reliability of the numerical results in this scheme. 

4.0.1 BLM summation 

Computing the leading-order 0{as) contribution does not tell us at which scale the cou- 
pling has to be evaluated; that is fixed only by terms of higher order in a^. The (extended) 
BLM prescription [35, 36, 37] is based on the conjecture that the bulk of higher order 
terms can be incorporated by replacing the fixed lowest order ccg by the scale-dependent 
as{k'^) where k denotes the momentum fiowing through the gluon line, cf. Fig. 1. The 
running of is determined by the first coefficient of the Gell-Mann-Low beta function: 
/3o = y^c — 1^/ — 9. The resulting expansion is of the type q;s(1 + Z^^i Cn(/3otts)") 
with n = 1,2,... representing the second-order, third- etc. BLM correction. The BLM 
prescription amounts to the assumption that the perturbative terms left out are much 

We show the distribution from Ref. [5] which did not incorporate 1 /ml subtraction corresponding to 
Darwin operator. 
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Figure 3: Gluon momentum scale distribution in rsi(6^c). Solid, dashed and dot-dashed lines 
correspond to /i = lGeV, 1.5 GeV and 2GeV, respectively; lighter short-dashed line illustrates 
the case of = (pole masses) . The area under each curve gives the first-order perturbative 
coefficient. 

smaller since not enhanced by the large value of /3o. A dedicated recent discussion can be 
found in Ref. [38]. 

This ansatz has turned out to be intriguingly successful in several cases. Often it is 
used in a simplified form by considering only the first nontrivial BLM correction - this 
constitutes the essence of the BLM scale fixing for [35]. However in applying this 
properly to semileptonic B decays one has to be aware of some specific complexities: 

• There are effectively two 'running' couplings in the problem, namely ag and nib- 

There is no reason why they should be controlled by the same scale - on the contrary, 
one expects them to be different. Moreover, the standard BLM scale fixing procedure 
is not defined - and rather often becomes meaningless - in the case of more than 
one coupling. Indeed, varying the normalization scale for masses the BLM scale for 
as varies in the whole range from zero to infinity. 

• As stated before good judgment has to be exercised in choosing the proper mass 
construction. The pole mass with its intrinsic infrared uncertainty ~ O(Aqcd) is 
ill-suited in this context. Using the BLM improvement only aggravates this problem. 

The analytic expressions for the perturbative corrections and the way to implement 
the Wilsonian cutoff are described in Appendices A. 2 and A. 3. They also give numerical 
values. 

4.0.2 Non-BLM contributions 

The non-BLM terms are expected to be small. They are known in the analytic form 
neglecting the charm mass (the case of Tsiib^u)) [39]. In 6— s>c a numerical evaluation 
exists. The authors of Ref. [40] have computed 0{al) corrections to 6 — > c£i/ at the three 
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different kinematical points = 0, and {nib — rric)'^, wfiere denotes tlie mass of 
tlie lepton pair. From tliese tliree values one interpolates over the full range. The non- 
BLM corrections indeed turn out to be small for the low-scale masses, which is a highly 
welcomed result: the consistent procedure to fix the scale of ag in the non-BLM terms is 
not known unless third-order corrections beyond simple BLM are available. This would 
yield significant numerical uncertainty if the non-BLM coefficients were not suppressed. 

4.1 Overall perturbative correction 

It is gratifying that the perturbative corrections to Tsi{b — >■ c) with the short-distance 
low-scale masses show good convergence [41]. Let us define 

zoir) A^'''\r; fx) = zo{r) + a^f + a2 {^f + as (^)' + ... (33) 

assuming the standard choice of the MS coupling ^^(M) normalized at M — mi,- As an 
example, at /i — l GeV, mft = 4.6GeV and mc = 1.15 GeV we have 

oi ~-0.94zo(r), 0^°"^-^^^ ~ -l.Ozo(r), 

af^^ ~ -0.45/3o^o(r), ^ ~ -0.21f3', Zoir), af ^ ~ +0.09/3o' ^o(r), ...(34) 

This can be compared with the expansion in terms of pole masses: 

ai ~-1.78zo(r), ar"^^"" ^ 1-4 ^o(r), (35) 
a^L^ ^ -1.92/3o zo{r), ^ ^ -2.8(3', zo{r), ^ ^ -A.9(3l zo(r), ... (36) 

(the second-order BLM coefficient in the pole mass scheme was first evaluated in Refs. [42] 
and its size interpreted as uncontrollable behavior of the perturbative series for the 
semileptonic widths). With this improvement, it fortunately does not make a too signif- 
icant difference which approximation beyond the two-loop result to adopt, as illustrated 
below. Using the complete BLM-resummed result has the advantage of being nearly in- 
dependent on the initial scale used to normalize ag. To this we add the second-order 
non-BLM correction evaluated with q;5 = 0.25 and regard this as our central estimate:^^ 

^P^'"*(r;/x) ~ 0.908 at = ^ and r = 0.0625 . (37) 

Alternatively, if we keep only the full second-order result evaluated with as{mb) = 0.22, 
this number becomes 0.909; including additionally third, or third and fourth BLM terms 
yields A^^^^ir; fx) ~ 0.903 and AP^'^ir; /i) ~ 0.904, respectively 

There are a few sources of possible theoretical uncertainties here. Before addressing 
them, we should emphasize that they shall not include uncertainties in the numerical 
values of as{mb) or in the running quark masses - those are determined from experiment 

^■^ Since the infrared domain is removed to all orders in perturbation theory, the effective scale cannot 
be too low. 
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and in principle have definite values as soon as the scheme is fixed, although they carry 
some error bars in practice. 

The uncertainty in the second-order non-BLM coefficient even estimated conserva- 
tively [40] leads to SA^^^^ about 0.0025; we may think it was actually overestimated. A 
potentially more significant uncertainty could come from varying the effective scale in as 
used to evaluate this term. However, since the magnitude of the non-BLM coefficient in 
Eq. (34) is about unity, this still would not exceed 0.003. 

The real question is then the size of the third- and even higher-order perturbative 
corrections. They are expected to be dominated by the ELM corrections which has 
been completely accounted for, and yield a sub-% correction. Hence, we have reasons 
to consider the estimate of the theoretical uncertainty in knowledge of the perturbative 
factor 5/lP'^'"*(r; 1 GeV) ^0.009 as conservative and justified: 

^pert^Q 252. 1 Qgv) = 0.908 ± 0.009 at Q;,(mb) = 0.22 . (38) 

The dependence on the actual value of as{mb) follows from Eq. (6). The dependence on 
the values of heavy quark masses and the used normalization scale /i can be obtained by 
explicit evaluation of equations in Appendices A. 2 and A. 3. 



4.2 Wilson coefficients of power-suppressed operators 

While the Wilson coefficient of the kinetic operator is identical to the coefficient of the 
unit operator ^p'^''* and therefore known completely to two loops, the ag corrections to 
the chromomagnetic operator have not been calculated so far. Once known, they would 
tell us the corrections to the coefficient for the LS operator as well. The coefficient for the 
Darwin operator is an independent one, and its perturbative corrections are not known 
either. Part of the perturbative corrections comes from the short-distance renormalization 
of the heavy quark vector and axial current, and this one always enters as an overall factor. 
Therefore, it is reasonable to factor out A^^"^^ in the expression for the width in the absence 
of explicit calculations of the remaining corrections. 

The accuracy of the tree level value of the exact chromomagnetic coefficient 1 + Ag^"^* 
is conservatively estimated as 30%. It is generally expected to deteriorate for higher 
operators, and we assume only a 50% accuracy for the Darwin coefficient; since the latter 
is still poorly known and is rather constrained from above by data [6], this does not 
introduce a significant additional uncertainty. The chromomagnetic expectation value, 
on the contrary is known rather accurately, Eq. (17). The actual normalization scale 
for total semileptonic width is probably somewhat higher, which slightly decreases the 
expected effective value of Hq in Eq. (1). We therefore estimate 

^ 0.01 , ^ ^ 0.015 . (39) 
Tsi Tsi 

Lack of evaluation of the (9(«s) corrections to chromomagnetic and Darwin Wilson coef- 
ficients becomes now one of the limiting factors in the overall theoretical precision [26]. 
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5 Theoretical Uncertainties in T(B — > Xcii^) 
5.1 First yet least: local duality violation 

A separate, more conceptual factor limiting the achievable theoretical accuracy would be 
violation of local duality. While it is related to high orders in the OPE, it has a different 
origin from what was discussed in Sect. 3.4 - a factorial growth of the expectation values 
of the higher-dimensional operators in the OPE, which makes the practical OPE series 
only asymptotic [27, 28]. Both the qualitative pattern and numerical aspects of duality 
violation depend crucially on the observable in question. 

In Ref. [43] a detailed and comprehensive analysis of limitations to local duality in 
B^Xciv was given. Here we reiterate some of the salient points: 

• Meaningful statements about duality violation require an accurate definition of the 
notion. That has been given for the OPE-treatable quantities hke inclusive decay 
widths, and must be understood as deviations from the nonperturbatively corrected 
predictions through a sufficiently high order in the OPE (as opposed to parton-level 
estimates) . 

• Claims of duality violations one can find in the literature often turn out to constitute 
violations of the OPE. No example of the latter has been given so far, and actually 
are not expected to exist. 

• On quite general grounds one infers that duality violations are given by 'oscillating' 
functions a Id sin that get further suppressed by powers of the energy release. 
No systematic excess or depletion of the rate compared to the OPE is allowed. 

• The fully integrated width rsi(S) relegates local duahty violation per se to rather 
high orders. I.e., limitations of local duality are suppressed by a high power of the 
energy release. However, a much larger impact can exist on differential distributions 
or various widths with cuts. 

• It has been estimated that limitations to duality in rsi(-B) amount to at most a 
few pcrmill. This bound cannot be rigorously proved from QCD without some 
quantitative understanding of details of dynamics governing the nonperturbative 
regime of QCD. However, in any scenario effects of local duahty violation in total 
semileptonic widths are much smaller than those in the differential distributions 
or moments used to extract the heavy quark parameters, thus introducing only 
negligible corrections in practice compared to the latter. 



5.2 Summary 

We have expressed the total semileptonic width of B mesons as a functions of heavy quark 
parameters, namely quark masses and B expectation values of local operators through 
order I/ttiq, in addition to The relation still has a potential theoretical error 
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due to pcrturbative uncertainties in the Wilson coefficient of tlie leading bb operator, 
the QCD radiative corrections in the Wilson coefficients of chromomagnetic and Darwin 
operators, in the impact of 'intrinsic charm' and in the contributions from higher-order 
power suppressed contributions. As already stated in Eq.(7), we estimate their impact 
on \Vcb\ not to exceed ±0.005, ±0.012, ±0.007 and 0.004, respectively. We think that 
these estimates are far from aggressive; the actual uncertainties might even be smaller. 
For instance, as mentioned in Sect. 6.1, there are correlations between and the value 
extracted for uiq. Furthermore, a large part of the contribution from the charm operators 
is accounted for by the Darwin contribution, when the latter is extracted from experiment, 
cf. Eq. (27). By the way, this "IC" contribution may explain why the central value of 
as suggested by experiment is smaller than anticipated. 

In summary: we estimate that the present theoretical accuracy in expressing rsi(i?) 
through the heavy quark parameters corresponds to an uncertainty in \Vcb\ of about 2%, 
Eq. (7). 

6 Determining the Heavy Quark Parameters 

Through order l/rriQ we have the following heavy quark parameters: 

1. Quark masses mbdi) and mdn); 

2. The chromomagnetic and kinetic expectation values Hc{n) and /^^(a*), respectively; 

3. The Darwin and LS terms p%{p) and Pls{i^)j respectively. 

We have excluded the charm operators represented by which are treated separately. 
We want to determine the values of these HQ parameters in a way that does not jeopardize 
the strong points of the OPE expression for the width, Eq. (1), namely an expansion in 
inverse powers of (at worst) m^ — rric, when only local operators are relevant. 

The b quark mass has been extracted from beauty production at threshold in e'^e~ 
annihilation by several authors [44]. Their findings are completely consistent within the 
stated uncertainties of about 1.5%. The techniques (and moments) employed in the 
analysis differ somewhat from author to author and the agreement in their findings is 
reassuring. One should keep in mind, though, that these determinations not only share 
their experimental input. The underlying approach is the same, as well as a number 
of general assumptions. The value stated in Eq. (15) could thus be subject to some 
systematic bias. Arguments based on the SV sum rules actually suggest that m^'"(l GeV) 
lies a bit above 4.6 GeV, if a relatively low kinetic expectation value extracted from 
experiment so far is confirmed. 

The values for nib and rric can be determined also from the shape of the energy spectra 
in semileptonic (or 6— >s + 7) B decays, which is most concisely encoded in their first few 
moments. Extractions were attempted based on the formulae of Rcf. [45] , their status is 
still controversial. Another relevant for us area of concern in the CLEO analysis is that 
the relation of Eq. (16) is imposed as a constraint for mb — rric in semileptonic decays. 
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Irrespective of that it should be noted that the central value extracted from e^e~ bb 
and from B decays agree well even within only the stated errors reported in the latter. 
Considering that the two processes have very different systematics experimentally as well 
as theoretically and should not have a correlation with the possible bias alluded to above, 
this is a nontrivial and encouraging result. 

Since the early days of the HQE one has often employed the expansion given in Eq. (16) 
that relates the difference mh—mc and the spin averaged charm and beauty meson masses. 
Yet in view of the accuracy in nib that has been achieved, the relation (16) may represents 
the most vulnerable part of the analysis. It suffers from two systematic weaknesses: 

• It brings in effectively an expansion in l/rric rather than l/m^ or l/(m;, — mg). At 
the very least it is obvious that l/(mft — mc)^ and 1/m^ differ a lot in magnitude, 
since a short-distance charm mass enters. 

• It involves the nonlocal correlators p^^ and intrinsically unrelated to B de- 
cays. There are actually complementing indications that their sizes are particularly 
large. ■'^^ 

Accordingly one is ill-advised to impose relation (16) in the analysis. Instead one should 
check the validity of the relation if there is an option of determining rric independently. It 
is possible that Eq. (16) holds with reasonable accuracy; for eigenvalues of Hamiltonian 
often exhibit a more obedient behavior than formfactors etc. However we cannot count 
much on such good behavior a priori. Since /i^ — /Xq-Ca*^, the expansion 

m,-m, = MB-Mn + (—-—)- 

2 Krric rrib/ 

(40) 

i.e., where the pseudoscalar rather than the spin averaged meson masses are used, is 
expected to be more stable with respect to higher order power corrections [14]. It requires 
some care incorporating the perturbative corrections to /jLq, but is more promising in 
reducing the uncontrollable errors. 

When the primary goal is to determine \ Vcb\ as precisely as possible rather than the HQ 
parameters, another observation becomes important. Almost the same combination of the 
HQ parameters controls the low semileptonic moments and the total semileptonic width. 
It also controls the first hadronic moment (Mx)- This combination of HQ parameters 
can thus be determined with higher accuracy than the individual HQ parameters - and 
the error associated with it becomes partially experimental. The theoretical status of 
evaluations of the moments is not as advanced as for Fgi, therefore the expressions used 
for them should be viewed with lower confidence than what we know about the width 
itself. Nevertheless, we shall adopt them literally to illustrate the point. 

^''it was independently observed analyzing the actual hadron mass spectrum [14], in the exactly solvable 
't Hooft model [46] and in the pilot lattice study [18]. 
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6.1 Numerical examples 

An important new element is provided by the measurements of Icpton energy and hadronic 
mass moments in inclusive semileptonic B decays, as pioneered by CLEO and achieved also 
by DELPHI. Such data provide us with novel detailed information on the HQ parameters. 
Even though they do not allow to pinpoint the values of and rric separately, they turn 
out to improve significantly the accuracy and credibility of extracting \Vcb\- To illustrate 
this, we adopt evaluation Eqs. (5) and (6) neglecting possible uncertainty in ^^(mb). 

Comparing with the predictions for, say, the first lepton moment (Ei) or first hadronic 
moment (Mj.) expressed in terms of the very same heavy quark parameters, we arrive at 

\Vcb\ = Kb {l - l.TOGeV"^ ((£;^)-1.383GeV) -0.075(mc(lGeV)-1.15GeV) 

+ 0.07 (/x^-0.4GeV2) - 0.055 (pf)-0.2GeV^) 

- 0.085 (/x|-0.35 GeV^) - 0.005 (^1^+0.15 GeV^) } , (41) 

\Vcb\ = V;°6 {l + 0.14 GeV-2((M|-) -4.54 GeV^) -0.03 (mc(lGeV)- 1.15 GeV) 

+0.1 (/i^-0.4GeV2) +0.1 (p^-0.2GeV^) 

- 0.01 (/x^-0.35 GeV^) + 0.006 (^pls+0.15 GeV^) } . (42) 

Taking for orientation the literal error bars on the two moments quoted by DELPHI [6] 
we would arrive at the uncertainties due to heavy quark parameters at the level of 2.5% 
and 1.2%, respectively. To have an independent experimental constraints on kinetic and 
Darwin expectation values one can use the second and third hadronic moments, or their 
improved version [26]. 

The additional uncertainty would come from the charm expectation value H^. How- 
ever, as stated by Eq. (27), in such an approach it is expected to be milder: extracting 
in practice the Darwin expectation value, one partially accounts in it for the effects of 
possible contribution of He- 

A similar reduction of sensitivity to a.^ looks probable. According to Fig. 1 the effective 
scale of heavy quark masses shaping rsi(6— >^c) is somewhat larger than 1 GeV; if the same 
applies to the lowest moments, there should be a correlation between the extracted values 
of the masses and as reducing the dependence on the latter. This can be visualized by 
adopting a higher normalization scale for masses //~ 1.5 or even 2 GeV which reduces the 
overall perturbative correction. 

7 Summary and Outlook 

The heavy quark expansion allows to express r(i? — > Xcii') as a series in powers of 
i / (nT-b—nT'c) with coefficients that are controlled by the B expectation values of local HQ 
operators; the number of unknowns here through order is limited. Duality violations 
for this fully integrated rate have been estimated not to exceed the few permill level, i.e. 
to be irrelevant. Contributions from unknown higher order corrections - perturbative as 
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well as nonpcrturbative - are also estimated to be near or below 1%. It would require a 
very major effort to improve on this situation; at present we see no need for that, except 
computing the perturbative corrections to the chromomagnetic and Darwin operators. 

The burning issue is how accurately we can determine the values of the HQ parameters 
from other observables. In the past the validity of the heavy quark expansion for r(S— >■ 
Xciiy) has been compromised by invoking the relation of Eq. (16): mf,— is effectively 
given by an expansion in l/rric with the additional drawback of uncontrolled non-local 
correlators contributing. 

This weak link in the analysis can be overcome now: measuring the energy and mass 
moments in B decays allows us to obtain values for the HQ parameters from expres- 
sions in powers of l/(m5 — m^), where contributions from nonlocal operators are absent. 
Furthermore the HQ parameters the size of which previously could be inferred mainly 
by theoretical arguments, can be measured - and actually in more than one way. Such 
redundant extractions are a powerful tool in establishing theoretical control. It turns out 
that the low leptonic and first hadronic moments are controlled by almost the same com- 
bination of HQ parameters as r{B ^ Xciiy) . This means that uncertainties previously 
viewed - correctly - as theoretical, become more experimental and thus can be reduced 
by better data. The situation there - while promising and encouraging - is still fluid re- 
quiring closer experimental as well as theoretical scrutiny. On the theoretical side higher 
order contributions to the moments have to be evaluated with the same care as it has 
been done in this paper for the total semileptonic B width. On the experimental side it 
is essential that the cut in the lepton energy be kept as low as possible. For unless the 
major part of the spectrum is measured, limitations to local duality can enter through the 
'back door' to haunt us. Furthermore such cuts decrease the hardness of the transition 
and thus deteriorate the convergence and reliability of the expansion. 

Extracting |Vc6| with neither the theoretical nor the experimental uncertainty exceed- 
ing the one percent level is thus within our reach. It appears unlikely that a higher 
accuracy should ever become necessary. 

Finally - and this is one of the strengths of the OPE - the same values of the heavy 
quark parameters m^, /x^, pi, etc. can be used when extracting \Vub\ from B 
transitions. The latter, however additionally require evaluation of the properly defined 
four-quark expectation values with h and light quarks. 
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Appendices 



A.l Heavy quark masses 

Having stated in Sect. 3.1 our preferences for which quark masses to adopt, we mention 
here how one can translate between the different 'quark mass languages'. 

The usually used MS mass mgdj,) is not really a quantity in the Lagrangian, it is 
rather a combination of parameters optimized for perturbative calculations in dimensional 
regularization. Its perturbative relation to the pole mass 



m 



pole _ TrrrfTrT 



m{m) s 1 + 



TT 



+ 



TT 
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is known now to three loops [47]. At n > mq the MS mass coincides more or less with 
the running Lagrangian mass probed at scale oc and thus is appropriately used to 
describe high energy processes, E ^ mq. However maintaining by definition the same 
//-dependence it becomes unphysical for // ^ mq, 



mqiii) ~ mq{mq) 



1 + 



2a, 



TT 



1 "^<9 

log— 



(A.2) 



where actual running becomes much slower. The MS mass is thus inappropriate for 
treating decay processes where the relevant scales are well below mq and the evolution 
to a low normalization point fi <C mq becomes crucial. 

The kinetic mass in perturbation theory order by order is defined by subtracting from 
the pole mass the perturbative contributions of the heavy quark parameters: 

pert 2mi){fl) 

where the latter arc determined from the corresponding SV sum rules cut at energy e = //: 

^2 13^ 
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with the standard notation A = limmQ-^[MHQ—mq], and 
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(A.5) 



where Cf = ^ and Ca = Nc = 3. and M denoting an arbitrary scale used to normalize as- 
The latter is assumed in the MS scheme in Eqs. (A.4) and (A.5). The thus defined mass 
exhibits only a linear in fi dependence on the scale, 

dm^^''(/x) _ 16 a,(/i) _ 4 a,(/i) /i 
d/i 9 TT 3 TT mq 



+ O {al) , 



(A.6) 
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compared to Eq. (A. 2). The pole mass would corresponds to integrating this evolution 
down to /i = at each particular order - alas it leads to an ill-defined result. 

One could have extended the definition shown in Eq. (A. 3) to include the corrections 
scaling like Oi^^ and higher. Besides the known perturbative contribution in this 

Q 

requires the similar terms for the nonlocal expectation values p^^ and p|; the second- 
order non-BLM pieces there have not been calculated. The one-loop effects, on the other 
hand, appear to cancel in the sum of all h^/itiq terms, and the resulting shift beyond one 
loop (say, including evaluated ELM corrections) is tiny. In any case for the beauty mass 
this makes no visible difference whatsoever. There is an additional more subtle reason not 
to pay much attention to such terms. Therefore, we prefer to consider Eq. (A. 3) as exact 
defining what is understood by the running masses, unless accuracy in their practical 
extraction will improve in the future by orders of magnitude. 

Often "HQET parameters" like A, — Ai etc. are quoted rather than A(/x), p^{p) with 
appropriate /x; they simply correspond to taking the limit /i — 0, and in this respect 
are the close relatives of the pole masses. Likewise they are not well-defined; moreover, 
the irreducible uncertainty and order-by-order instability in their values is of order unity 
compared to their magnitude. In concrete calculations it is often possible to identify their 
relation to well-defined Wilsonian expectation values. Yet one has to be aware that their 
values depend sensitively on the context in which they appear, and one cannot transfer 
them from one case to the next blindly. As a rule of thumb one has for routinely used 
first-order calculations with the second-order ELM improvement using as{m}}) coupling 
equated to 0.22 



Ahqet ^ A(l GeV) - 0.255 GeV , 



-Ai ~ pI{1 GeV) - 0.18( GeV)^ . 



(A.7) 



Finally we quote the relation between the running mQ{p) and MS mass mQ{mQ) 
including full and third-order ELM corrections: 



\ , 4 a^S(m,) 

3 TT 





3 rrib 



71 



2m2 



^ I ^ + ^ 1 +^ + ^(21n2-^ 



48 



, rrih H 

m h- 

2fi 3 




8^ 

9 \ 2fi 6/ y \) I mt 

3(2353 13 , 7 16 ^ 

[2592 36 6^^ ' 9 




mb 13 

in 

2/i 6 



8 




67 



TT 



In— H — H 

2fi 37 36 6 



Ik. 

3 w? 



In 



ruh 
2ii 



vr 



(A.8) 



This can serve as a useful reference to translate between the kinetic and other possible 
low-scale running masses, along with the similar relation between the pole and kinetic 
masses. 
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A. 2 Perturbative corrections in the Wilsonian approach 



Here we give relations allowing one to determine in practice perturbative coefficients in 
the approach employing the Wilsonian separation of large and small virtualities. 

The starting point is the general short-distance expansion exemplified by Eq. (8) for 
^si{b c), or its practically employed form Eq. (1). As stressed in Sect. 2, it must 
be //-independent. Then one considers it in perturbation theory limited in practice to 
a particular order in and l/mt. Requiring yU-independence completely defines the 
//-dependence of the Wilson coefficients once the /x-dependcncc of masses and matrix 
elements in perturbation theory is known [9, 12]. For truncated perturbation theory the 
actual expression has a residual dependence on /i to the uncalculated orders in ctg. On the 
other hand, in the finite-order perturbation theory one can formally put /i — O. This would 
set all the power-suppressed expectation values to zero and yield pole masses, hence fixing 
the perturbative series for A^'^^^{r; 0) to be the 'conventional' perturbative corrections in 
the pole mass scheme. This provides the initial condition for perturbative /i-dependence 
at fjL — O. 

To give the practical implementation for rsi(6— >c) we need to specify more precisely 
our definition in Eq. (33): now we write 

zo{r) ^P^'-*(r; //) = zo{r) + ai(r; //)f + a^ir; //) )' + as{r; //) )' + ... (A.9) 



and use the short-hand notations a 
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fjL — O (pole mass scheme) . Then one has 
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and the coefficient 0^°^ (r) is given by [48, 49] 

af\r) = ^ [75 - 12 77^- 956 r- 192 7rV^ + 956 r^- 3 (25 + 4 TT^) r^ 
+ 3847rVt (l + r)+4 (-17 + 64r - 64r^ + 17r^) In(l-r) 
-4r (60 + 270r-4r^ + 17r^- 384Vr (1 + r) ln(l + ^/f)) ln(r) 
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+ 48 (l-16ri + 30r2-16ri + r^) In(l-r) lnr-12r^ (36 + r^) InV 

3 / 128 r5 fl + 
- 3072r2 (1 + r) Lial^r ) + 18 4 + 64r^ + 4r^ + ^ '- \ Li2(r) 



The combinations Ai 2, Pi,2 and pi_2 are the coefficients of the perturbative expansion of 
[^(/^)]pert and [/x^(/xj]^ert in Eqs. (A.4) and (A.5) and of [pUfJ')]^erv ^ first and 



second orders in as, respectively: 
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It should be noted that we did not include here shifts due to the spin-singlet l/mt 
pieces of the kinetic and chromomagnetic operators (we have computed them to order as 
and to all orders in BLM). These corrections are 1/mf and governed by the full l/2mb 
scale, hence totally insignificant in practice. We did include the effect of the Darwin 
operator to two full loops (modulo uncalculated 0{as) piece A'j^^^), since its coefficient is 
enhanced and strongly dominates over the effects. 

The improvement by the Wilsoiiian cutoff of the resulting perturbative series can be 
illustrated by comparing them ai p — 1 GeV 
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with those in the 'pole' (HQET) scheme: 
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where numbers refer to mc/m^ = 0.25. 



A. 3 BLM summation with Wilsonian cutoff 



The BLM correction in practice are simpler to calculate directly, using the generahzed 
order-Q;^ relation (A. 10) between Wilson coefficients with and without the cutoff. The 
technique of BLM resummation was recently reviewed in this context in Ref. [5]. It 
requires computing one-loop corrections with a non-zero gluon mass A. The expressions 
for the one-loop correction to the perturbative width with non-zero gluon mass are rather 



lengthy and can be found in Ref. 



One also needs similar terms in 
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and [pf)]'''^'^*- These are calculated by integrating the SV spectral density of Ref. [12] with 
the proper power of energy and are given as follows: 



Ai(/x; A) 
One then has 
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The perturbative coefficients are exphcitly given by the series 
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and the resummed result reads as 
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Numerically at mc/mb — 0.25, ///m6 = l/4.6 and as{mb) — 0.22 we obtain 
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